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Abstract 

We analyzed the effects of population density and clinnatic variables on the rate 
of population growth in the guanaco (Lama guanicoe), a wild cannelid species in 
South Annerica. We used a tinne series of 36 years (1977-2012) of population 
sannpling in Tierra del Fuego, Chile. Individuals were grouped in three 
age-classes: newborns, juveniles, and adults; for each year a fennale 
population transition nnatrix was constructed, and the population growth rate (A) 
was estinnated for each year as the nnatrix highest positive eigenvalue. We 
applied a regression analysis with finite population growth rate (A) as 
dependent variable, and total guanaco population, sheep population, annual 
nnean precipitation, and winter nnean tennperature as independent variables, 
with and without tinne lags. The effect of population size was statistically 
significant, but the effects of the sheep population and the clinnatic variables on 
guanaco population growth rate were not statistically significant. 
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Amendments from Version 1 

We evaluated the comments of both referees (Dr. Baldi and 
Dr. Acebes) and considered them to be very useful. Following 
the reviewers' suggestions we: 

1) Improved the description and justification of the methodology 
used with respect to the guanaco sampling. In relation to the 
estimation of the guanaco population size we specified the 
use of the method proposed by Leopold, and we developed 
the necessary assumptions for this method to be valid. Brief 
comments about the surveyed area were also included to 
support the sampling methodology. Some questions related 
to the population sampling with the Distance method were not 
considered because we didn't use that method. 

2) Re-processed the regressions incorporating the sheep 
population as a new independent variable of the guanaco 
population growth rate {X); so our new analysis now has 

6 independent variables: sheep population, guanaco 
population size, winter temperature, winter temperature with 
a lag of 1 year, annual precipitation and annual precipitation 
with a lag of 4 years. Our main results didn't change and it 
was confirmed that only the guanaco population size seems 
to have an effect on the population growth rate {X). 

3) Complied with Dr. Acebes suggestion to use the Akaike 
measure of goodness of fit and so we changed the software 
used in the analysis (from Statistica to language R). 

As a consequence of these changes in the analysis. Table 1 and 
Table S1 (now referred to as Data set 1) were modified. 
Lastly, owing to these changes and corrections some citation were 
removed and others were incorporated. 

See referee reports 



Introduction 

In order to understand population dynamics and optimize the man- 
agement of wildlife populations it is important to identify how 
populations are affected by environmental factors and their own 
density, particularly in species living in extreme environments. 
Mechanisms that regulate populations of the guanacos (Lama 
guanicoe) are poorly known; this species was heavily exploited in 
the Patagonian steppe, and its numbers and distribution have dimin- 
ished significantly during the last century because of grazing con- 
flicts with a sheep-based society, and overhunting (the guanaco 's 
historical distribution has been reduced by 75% in Chile and Peru, 
and by 60% in Argentina) ^ The guanaco is in the Least Concern 
lUCN Red List category (http://www.iucnredlist.org/, downloaded 
on 11 September 2013) but it has been included in Appendix II 
of CITES 2013 (http://www.cites.org/eng/app/2013/E-Appendi- 
ces-2013-06-12.pdf). 

No population of any species can grow indefinitely, and population 
checks based upon different processes restrict population size and/ 
or geographic distribution; these processes are either density- 
stabilizing or density-limiting-. The former are of a biotic nature 
and depend on the interaction between individuals of the same or 
different species, while the latter are independent of population 
size. Stabilization results from density-dependence, with a regula- 
tory effect that varies in intensity with the size or density of the 
population itself, however, not all density-dependent factors are 
density-stabilizing. The density-limiting factors can also be consid- 
ered density-responsive because the per capita amount or availabil- 
ity of resources decreases as the population density increases. Thus, 
these two types of factors (density- stabilizing and density-limiting) 



rarely act independently: the density-limiting factors (generally 
of a physical and/or climatic nature), may determine the level at 
which populations become stabilized by the density-dependent pro- 
cesses, but they do not have a stabilizing capacity per se\ for this 
reason many wildlife population dynamic and management mod- 
els include the effects of climatic variables (e.g., Dennis & Otten, 
2000; Colchero et al, 2009)''l 

The guanaco is one of the two extant wild South American came- 
lids, and ranges from Northern Peru through Chile, and across 
Patagonia to southern Argentina and Chile, reaching Tierra del 
Fuego. It is found from sea level to nearly 4500 m on the Andes 
mountain range, and occupies a wide variety of habitats from hard- 
pan deserts to scrublands to grasslands^ 

Although there are several studies on the population regulation pro- 
cesses in mammals in generaF''^, and in ungulates in particular^' 
there are very few studies in wild South American camelids that 
analyze population growth regulation. There are several studies 
dealing with population dynamics of vicuna (Vicugna vicugnay^'^"^ 
and guanacos^^"^^; however none of them considered the effect of 
density, livestock and climate on population growth rate. 

Thus, our purpose in this work was to test the hypothesis that popu- 
lation density, sheep stock, and climatic variables, affect the growth 
rate of the guanaco population from the island of Tierra del Fuego, 
Chile. 

Materials and methods 

Study area 

We analyzed the guanaco population of the "Cameron" ranch (-53.9 S, 
-69.3 W), with an area of the 2000 km^ located in the Southern 
region of the Tierra del Fuego island, Chile. The altitude range is 
0-300 m above sea level, and the ranch is a mosaic of steppe with 
fragments of forest; the latter is a deciduous forest (Nothofagus spp) 
while the steppe is composed of meadows, peats and prairies. The 
forested area is about 8.8% of the total study area, with some degree 
of forest clearance that offers adequate visibility for guanaco sam- 
pling. The guanaco is the dominant wild herbivore, while sheep are 
the dominant domestic species; densities of the latter have fluctu- 
ated in the last decades between 1 1 and 23 sheep/km^^l In contrast 
to the mainland, the puma (Puma concolor), the main predator of 
guanacos, is absent on the island; since 1977 guanaco hunting has 
been controlled by the Chilean National Forest Corporation (CONAF, 
according to its Spanish acronym) and the Chilean Agricultural and 
Livestock Service (SAG, according to its Spanish acronym). 

The climate of the Tierra del Fuego island is characterized by an 
average annual precipitation of 200-400 mm, while in the forested 
areas the average precipitation fluctuates between 400-600 mm per 
year^'^. The average annual temperature is 6.5 °C and the average 
winter temperature is 2.2°C. 

Guanaco population sampling 

Guanacos were counted for 34 consecutive years between 1977 
and 2012 (except in 1986 and 1996), using the transect method 
with an undefined band width, from 1977 to 1995, and with a fixed 
width band from 1996 to 2012 (the latter with a maximum of 500 m 
to each side of the transect)^" --. The sampling methodology was 
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changed by expert recommendation made in 1995^^; in 2000 the 
fixed and the undefined band width methods were appUed simul- 
taneously, we tested both results and found no significant statisti- 
cal difference between them. The sampling period was carried out 
in the autumn and lasted approximately 7 days between 10:30 and 
19:00 h, with two observers in two 4x4 vehicles going over the 
main, secondary and local roads at speeds slower than 40 km/h. 
Despite randomly selected transects being recommended^^, pre- 
existing roads were used because according to Soto (2010)^^ the 
existing system of roads is an adequate sample of the area. Each 
road was covered only once, and in addition to individual guanaco 
counts, the following were recorded: weather conditions, time, dis- 
tance (km) from the starting point, GPS coordinates, observation 
distance from the transect (m) to guanacos, an estimate of the angle 
to the animal's position, and - when the animals were observed in 
groups - the number of individuals, the type of social group, and 
its structure (sex, and the age class: newborn, juvenile and adult). 
Observation distances were made by naked eye. The road network 
and all geo-referenced observations were processed with the Arc 
View 9.3 Geographical Information System (SIG), and transferred 
using program Map Source®. The cartography was kindly provided 
by the Chilean Agricultural and Livestock Service (SAG). 

The population size was estimated as given in Soto^^ which was 
based on the King method modified by Leopold (1933)^^ described 
by Raedeke (1978)^^ given by: 



N = 



A ^ n 

2 ^ y * X 



(1) 



where N is the total population to be estimated, A is the total study 
area surface, n is the total number of animals counted, x is the total 
transect distance covered (m) rounded to one meter, and y is the 
average of the perpendicular distance (m) from the transect to the 
animals counted (the factor of 2 is included to consider that there 
is one band to each side of the transect). The variance (S^) was esti- 
mated by (2), with p = n/N, and used to estimate the 95% upper and 
lower confidence intervals. 



p2 



1-p + n 
n + 2 



(2) 



Raedeke (1978)^^^ claims that the Leopold method is valid when the 
following conditions are satisfied: (1) the road systems must be an 
adequate sampling of the study area, (2) the network road must be 
randomly distributed, and (3) the animals included in the sampling 
must be randomly distributed in relation to the observer's route and 
he considers that in the south of the Tierra del Fuego Island these 
assumptions are fulfilled. Soto 2010^^ compared population estima- 
tions by the Leopold and Distance methods and found that the value 
of the means estimate by Leopold methods fall within the confi- 
dence intervals estimated by Distance; additionally, as all sampling 
periods used the same methodology whatever bias might exist in 
the estimation of the mean abundance using the Leopold method, 
the relative changes among years (and thus the temporal trend) will 
be adequately represented; thus, we conclude that the Leopold's 
method seems adequate for our purposes 



The area effectively surveyed in each sampling period was around 
420 km^ (about 20%, of the area under study). 

Climatic variables 

We evaluated the possible impact of two climatic variables on gua- 
naco population growth rate: average annual precipitation (mm/ 
year) and winter temperature (as the average of the temperatures of 
June, July and August, in °C). These climatic variables were selected 
because some studies suggested that they have an influence on 
guanaco demographic parameters: Rey et al. 2012-^' observed that 
after a severe drought the proportion of guanaco yearlings/females 
decreased significantly, and Sarno et al. 1999^^ detected a negative 
effect of winter snowfall on guanaco juvenile survival; because we 
had insufficient snowfall data we used the winter temperature as a 
proxy for winter snowfall. We used the 25 years (1977-2002) pre- 
cipitation and temperature time series of the CRU TS 2. 1 database, 
compiled by the Tyndall Centre, Climatic Research Unit, School of 
Environmental Sciences of the University of East Anglia, United 
Kingdom (http://www.cru.uea.ac.uk/cru/data/hrg.htm). As the CRU 
TS 2.1 data ended in 2002, we completed that time series for 2003 
to 2012 from the closest meteorological station to the Cameron 
ranch: Punta Arenas (Chile); this data was downloaded from the 
Internet site of the Meteorological Service of Chile (http://www. 
meteochile.gob.cl/). 

Sheep population 

The livestock effect was evaluated by considering the sheep popula- 
tion, based on time series data obtained from Soto (personal com- 
munication). Only eight years of data were available (1980, 1985, 
1990, 1995, 2000, 2005, 2008 and 2011); as the sheep population 
change between years was relatively smooth, we interpolated lin- 
early between two consecutive data points to generate a complete 
time series of 36 years. 

Guanaco population analysis 

We estimated the finite rate of increase (X) as a measure of pop- 
ulation growth rate for each of the 36 years of the guanaco data 
(population estimates were not available for years 1986 and 1996, 
and were linearly interpolated). To estimate X we developed a three 
age-class (newborns, juveniles and adults) female-only matrix 
model, and estimated the matrix parameters by fitting the model 
predictions to the field total population estimates; we summed the 
three age classes and multiplied the result by two (guanacos have an 
approximately 50% sex-ratio). The fit was carried out in an Excel® 
spreadsheet using the Solver tool, and the sum of squares (SSQ) as 
goodness of fit criterion (details can be found in Rabinovich and 
Zubillaga, 2012^^). This process resulted in a set of 36 population 
age-class structured matrices, and from each matrix we calculated 
the largest positive eigenvalue-^ as an estimate of X, using the 
PopTools add-in, an Excel program developed by Greg Hood (http:// 
www. c se . c siro . au/poptools/) . 

Population analysis 

To test our hypothesis, we used a multiple regression analysis to 
check the relationship between the population growth rate (X) as a 
dependent variable and guanaco population size (in natural logs), 
sheep population size (in natural logs), and climatic variables (annual 
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mean precipitation in mm/year, and winter mean temperature in °C) 
as independent variables. In the case of the chmatic variables we 
also evaluated the effects of various time lags (T), with T = 1 to 
7 years for precipitation, and T = 1 year for winter temperature; 
the lags were applied by averaging the previous T lagged years, as 
suggested by Shaw et al. (2012)^^ for the vicuna population (these 
lagged variables were used in the regression analyses in addition 
to the non-lagged precipitation and winter temperature variables). 
The inclusion of time lags is convenient because it may detect 
possible cumulative climatic effects on ungulate population growth 
rates as has been observed in, for example, deer and moose^*^ and 
in vicunas ^\ The multiple regression was carried out using the sta- 
tistical package "/m" of the RStudio software (Version 0.97.449 - 
©2009-2012 RStudio, Inc.). Since for multiple regression analysis 
it is necessary to check that there is no coUinearity (significant 
correlations among the predictor variables), we carried out a cor- 
relation analysis between independent variables, and used a p value 
< 0.05 to determine statistical significance. We used the Akaike 
index (AIC) as a goodness of fit criterion. 

Results 

Data set 1 shows for each of the 36 years of guanaco sampling all 
the variables (dependent and independent) used in the regression 
analysis. 

Dataset 1. Guanaco population abundance, sheep population 
abundance, estimated population annual rate of growth, average 
annual precipitation, and average winter temperature of the 
Cameron ranch, Tierra del Fuego. 

Guanaco population abundance (individuals), lower and upper 
limits of the 95% confidence intervals (CI), sheep population 
abundance (individuals), estimated population annual rate of growth 
(lambda), average annual precipitation (P; mm/year), and average 
winter temperature (celsius). Guanaco population values were not 
available for years 1986 and 1996, and were linearly interpolated; 
sheep population values were only available for years 1980, 1985, 
1 990, 1 995, 2000, 2005, 2008 and 201 1 ; other values were linearly 
interpolated. 

Click here to access the data. 

http://dx.d0i.0rg/l 0.5256/f 1 000research.2375.d27483 

V J 

Correlation analysis between the climatic variables showed a statis- 
tically significant correlation {p < 0.05) only among the precipita- 
tion variables for all the 1 to 7 year lags (but not with non-lagged 
precipitation), while the correlation results between all other inde- 
pendent climatic variables were not statistically significant. Due to 
the statistically significant correlation among the lagged precipita- 
tions we carried out seven independent regression analyses, all of 
them with the same independent and dependent variables, except 
the independent variables of lagged-precipitation, which were used 
in each regression run with a different lag value. We observed that 
the regression with precipitation lag T = 3 was the only one, aside 
from total guanaco population size, that had the intercept statistically 
significant. In all case the multiple regression analysis indicated 
that the total population (L/iNtot) was the only statistically signifi- 
cant variable (AIC = -81.599 p < 0.0514 and F = 2.414) (Table 1). 



We also carried out a simple regression analysis using X as the 
dependent variable and only the guanaco population size (in natural 
logs, LnNtot) as an independent variable, in order to compare this 
result with the full model (guanaco population, sheep population 
and climatic variables as independent variables). Table 2 shows the 
result of this second regression analysis, which was highly signifi- 
cant (AIC = -87.847, p < 0.0018, and F = 1 1.48,) and the resulting 
regression equation was X = 1.6373 - 0.0552 (±0.0163) L^Ntot. 

The negative slope of the regression equation is an indication of a 
density-dependent process, since when the population increases the 
per capita population growth rate decreases (Figure 1). The intercept 
of the regression line at A = 1 (population at equilibrium) results in 
a population size of 103,000 guanacos (51 guanaco/km^), which 
can be considered as an estimate of the carrying capacity of the 
Cameron ranch for this guanaco population. 

Discussion 

We were not able to confirm the expected effect of either climatic 
variables or sheep population on the population rate of growth (X) 
in the guanaco population of the Cameron ranch. We expected an 
effect of climatic variables because Samo et al. (1999)^^ found a rela- 
tionship between winter weather and guanaco yearling's survival, 
and Rey et al. (2012)^^' recorded the effects of a drought on guanaco 



Table 1. Results of multiple regression analysis between the 
population growth rate (k) as dependent variable, and (in log 
scale) total guanaco population (LnNtot), sheep population 
(LnSheep), annual mean precipitation (Annual Precip), mean 
precipitation with lags of 3 years (Precip (T = 3)), winter mean 
temperature (Winter temp.), and winter mean temperature with 
1 year of lag (Winter temp (T = 1)) as independent variables. 
The results show the regression coefficients (e Estimate), their 
standard errors (Std. Error) and t-test values (t value), with their 
probability value (Pr(>ltl)). 





Estimate 


Std. Error 


t value 


Pr(>ltl) 


Intercept 


1.8840 


0.8682 


2.1700 


0.0383 


LnNtot 


-0.0492 


0.0220 


-2.2370 


0.0331 


LnSheep 


-0.0272 


0.0654 


-0.4150 


0.6809 


Annual Precip 


-0.0001 


0.0002 


-0.3580 


0.7228 


Winter temp 


0.0173 


0.0191 


0.9030 


0.3738 


Precip (T = 3) 


-0.0003 


0.0005 


-0.5590 


0.5803 


Winter temp (T = 1) 


0.0309 


0.0182 


1.6980 


0.1002 





Table 2. Results of the simple regression analysis between the 
population growth rate (k) as dependent variable, and (in log 
scale) total guanaco population (LnNtot). Interpretation of the 
statistical parameters as in Table 1 . 





Estimate 


Std. Error 


t value 


Pr(>ltl) 


Intercept 


1 .6373 


0.1624 


10.082 


9.45E-12 


LnNtot 


-0.0552 


0.0163 


-3.388 


0.00179 
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9.5 10.0 10.5 11.0 11.5 
LnNtot 



Figure 1 . Effect of the guanaco population size on guanaco's finite population growth rate {k). On the x-axis the guanaco population from 
the Cameron ranch (Tierra del Fuego, Chile) was transformed to a natural logarithmic scale. The values of X represent the finite population 
growth rate (on a per year time unit). 



fecundity. Thus we anticipated that as winter temperature and/or 
annual precipitation decrease, the guanaco population growth rate 
might also decrease, as shown by Mech et al 1987-'^ with deer and 
moose. They also used a simple linear regression, and found that 
as winter snow accumulation increased, annual percentage change 
in population numbers decreased. With respect to sheep popula- 
tion, we also expected to find some type of relationship because 
this domestic herbivore is considered a competitor of the guanaco^^^ 
for food and water. Our negative results are even more surprising 
because in the last 10-12 years of the time-series data the guanaco 
population fluctuated remarkably, suggesting some extraneous 
factor in addition to between population density, maybe climatic 
variables. It is well known that environmental effects are more 
important when the population size is near the carrying capacity, 
particularly in large mammals characterized by low reproductive 
rates, long life-spans, and populations that are resource-limited 
(features typical of species referred to as the ''K-selected" species)^ ^ 
The difference between our results and those of Samo et al (1999)^^ 
may reflect that these authors used the average winter snowfall 
while we considered average winter temperature. 

The results of the guanaco population sampling suggest a certain 
trend in the population size, with a more or less exponential growth 
in the first few years, becoming more variable as the population 
grows, with more marked fluctuations during the last 10-12 years; 
this may imply that the population is becoming stabilized, possibly 
approaching its carrying capacity. On the other hand, the sampling 
results of those last 10-12 years show abrupt "jumps" in some years 
(for example 2004-2005) that may be the consequence of the mobil- 
ity of guanacos from the Cameron ranch to neighboring sites and 
vice versa; since the transect sampling does not identify individuals 
the effects of possible local displacements could not be considered. 

On the other hand the lack of effects of the sheep population size 
on the finite population growth rate (X) in this guanaco population 



conforms well with a recent study^- in Southern Chile that found 
that the potential competition between guanaco and sheep is low. 

In contrast to our results, a study by Taper & Gogan (2002)^^ on 
the Yellowstone elk population, with weather variables included as 
covariates within a dynamic model, found that spring precipitation 
had a positive regression coefficient. 

Contrary to what was expected based on the literature of ungulate 
population dynamics, weather variables do not seem to influence 
the density-dependent population growth rate of the Cameron ranch 
guanaco population. Our aim was to carry out a preliminary anal- 
ysis that would help identify some of the regulatory processes in 
this guanaco population and so we used a simple multiple regres- 
sion analysis as a first step in that direction. However, we note that 
in order to test the effects of climatic variables on population regula- 
tion of large ungulates such as the guanaco, the use of a population 
dynamic model would be recommended. A population dynamic 
model would better account for the interaction between density- 
dependent processes and weather variables than a simple multiple 
regression between the latter and guanaco population growth rate. 

Data availability 

FIOOOResearch: Dataset 1. Guanaco population abundance, sheep 
population abundance, estimated population annual rate of growth, 
average annual precipitation, and average winter temperature of the 
Cameron ranch, Tierra del Fuego, 10.5256/f 1000research.2375.d27483 
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Referee Report: 21 August 2014 
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The authors have improved the ms from their first submitted version. Please see the comments below to 
further improve the article: 

• First paragraph (introduction): please put guanaco in the singular. 



• Guanaco population sampling section: please replace 'SIG' with 'GIS'. 

• Guanaco population sampling section: please remove 'surface' in "A is the total study area 
surface". 

• Results: part of the results is a justification of the analyses carried out and should be placed in 
Material and Methods. 

• Discussion: in my view one of the most important results is the absence of negative effects of 
sheep on guanaco population growth rate that indeed has been cited as one of the major drivers of 
guanaco population decline particularly in Patagonia in the 20 century. For that reason I think 
authors should discuss this important result more deeply. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 



Competing Interests: No competing interests were disclosed. 
Gary Luck 

Institute for Land, Water, and Society, Charles Sturt University, Albury, Australia 



Approved: 04 August 2014 

Referee Report: 04 August 2014 
doi:1 0.5256/f 1 0OOresearch.4501 .r2049 

I have cross checked the comments from the two original referees with the responses from the authors 
and consider that the authors have adequately addressed all comments. The study has various 
limitations, such as a limited set of predictor variables, lack of the use of a sampling technique (e.g. 
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Distance sampling) that accounts for differences in detectability, and application of a simple regression 
analysis. The authors have adequately addressed each of these major issues in their response to 
referees - though readers should be aware of these limitations when interpreting the results of the study. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 

Competing Interests: No competing interests were disclosed. 



Referee Responses for Version 1 

I gy I Pablo Acebes 

^ • Department of Ecology, Autonomous University of Madrid, Madrid, Spain 

Approved with reservations: 19 December 2013 

Referee Report: 19 December 2013 
doi:1 0.5256/f 1 000research.2588.r2847 

The authors explore the population dynamics of a long-term monitored guanaco population (36 years) in 
the Chilean Tierra del Fuego, trying to discern if the process is governed by population size and/or 
environmental (climatic) variables, assuming that guanaco growth rate follows a density-dependent 
process. Although the topic is very interesting, the authors have room to improve the manuscript. Below 
the authors will find major and minor comments. 

Major comments 

The manuscript is classified as a "Research article" when indeed it seems a "Short note". It would greatly 
benefit from deeper and more detailed explanations of: 

1 . The theory of ungulate population dynamics, driven by density-dependent and 
density-independent processes; please expand this topic in the introduction and discuss more 
deeply your results considering the theoretical framework (see comments below). 

2. The census methodology (see comments below). 

3. Statistical modeling (see comments below). 

Ungulate populations are regulated via endogenous forces (density-dependent processes), and 
exogenous forces (density-independent): weather and predation. Weather determines the primary 
production, i.e. the resource (forage) availability, which in turn determines the carrying capacity. On the 
other hand weather, through extreme conditions (e.g. severe winter climate), modulates population 
dynamics (starvation, survival). Predation is considered an important force which shapes population 
dynamics as well. In addition, poaching or overhunting act as predation, and competition with other 
ungulates is another important exogenous force. With this framework authors should justify: 

1 . The background supporting the study, as it currently is a bit poor, in my view. The authors should 
explain more in detail, for example, which are the differences between the density-stabilizing and 
the density-limiting processes, or the evidence of including climatic variables when modelling 
population growth rate trends. 
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2. The use of these climatic variables: the average annual precipitation and winter temperature. 
Without any justification, authors could have also included others variables in the models, as 
climatic variables do not account for the pattern described, but this doesn't necessarily mean that 
climate hasn't got any effect on the guanaco population growth rate. Moreover authors should 
explain if these two climatic variables are selected as indicators of primary production and harsh 
weather conditions. 

3. As authors stated in introduction, guanaco populations have undergone overhunting (poaching) 
and competition with sheep, these being two of the major forces of their dramatic decline in the last 
century. For that reason, it is quite surprising that no variables related to livestock were included in 
the modelling procedure. The population has increased a lot in the studied period (x8) and just 
climate and/or density-dependent process do not explain alone such increase in my view. 
Therefore authors should consider including livestock effects, its influence on resource availability 
along the studied period and the mentioned overhunting, as this could really explain the detected 
trend. 

The guanaco sampling method has changed during the studied period (different band widths). Besides, 
were the transects chosen based on the topography and vegetation cover of the ranch to ensure that all 
areas were visible to census-takers along their transects? I am not sure if that could be a source of error, 
but authors should discuss this point. Moreover, I recommend authors to re-analyze their data using the 
DISTANCE program (Buckland etal., 2004) to estimate population size over the period, as the analyses 
performed by this package seems much more robust that the one employed by authors. If they do so, 
authors could also include topographic variables in the models to see whether the resulting models 
change. The reason is that the authors do not explain if guanacos occur in both biomes (steppe and 
forest), because if that was true, detection of guanacos was completely different in both ecosystems. 
Therefore the band width of 500 or 1 000m is not valid for forests. 

Regarding the population size estimated along the 36 years (Supplementary table), important (huge) 
differences are found, for example between 2004 and 2005 or between 2006 and 2007, or between 2007 
and 2008. My concern is if such "jumps" could be related to the census surveyed. Please discuss this 
point. 

The analytical approach employed by the authors doesn't seem appropriate in my view, as they also 
recognize in the discussion, and that could be the reason for not finding any climatic variable significant in 
the model: 

1 . The authors use a stepwise regression when other approaches such as logistic regressions, Ricker 
or Beverton-Holt models could give more suitable results. Please explain or justify your decision. 

2. The selection of the best model among the candidates should be done by Akaike's information. 
Minor comments 

1 . Refer to climatic variables or climatic effects, instead of climatic covariables. 

2. Explain more in detail Corani & Gatto's studies. 

3. Please remove information related to metabolic studies between guanacos and sheep (1st 
paragraph of study area section). 

4. Please add Nothofagus spp among the brackets after deciduous forest (study area section) if it is 
the case. 
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Give more information about livestock rancliing, density, and size and the relation of livestock 
rangers with guanacos along the studied period (study area section). 



5. 



6. Are the guanacos found both in the steppe and the forest biomes? 

7. Please re-write the paragraph related to guanaco sampling through transect method in the 
consecutive periods (it is not clear enough). 

8. Please remove information related to binoculars and GPS brands (did the authors use the same 
gadgets across 36 years?). 

9. Explain the interest of using time lags (T) as covariables in the model. 

1 0. Please add the Statistica version employed. 

1 1 . Please move explanations of the analyses perform from results to material & methods section. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard, however I have significant reservations, as outlined 
above. 

Competing Interests: No competing interests were disclosed. 
1 Comment 

Author Response 

Maria Zubillaga, CONIGET-UNLP, Argentina 
Posted: 14 May 2014 

• "The background supporting the study, as it currently is a bit poor, in my view. The authors 
should explain more in detail, for example, which are the differences between the 
density-stabilizing and the density-limiting processes, or the evidence of including climatic 
variables when modelling population growth rate trends. " 

After the following paragraph that was already in the original text: 
"No population of any species can grow indefinitely, and population checks based upon different 
processes restrict population size and/or geographic distribution; these processes are either 
density-stabilizing or density-limiting" 

We added the following paragraph explaining the differences between the density-stabilizing and 
the density-limiting processes. 

"The former are of a biotic nature and depends on the interaction between individuals of the same 
or different species, while the latter are independent of population size. Stabilization results from 
density-dependence, with a regulatory effect that varies in intensity with the size or density of the 
population itself, however, not all density-dependent factors are density-stabilizing. The 
density-limiting factors can also be considered density-responsive because the per capita amount 
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or availability of resources decreases as the population density increases. Thus, these two types of 
factors (density-stabilizing and density-limiting) rarely act independently: the density-limiting 
factors (generally of a physical and/or climatic nature), may determine the level at which 
populations become stabilized by the density-dependent processes, but they do not have a 
stabilizing capacity per se; for this reason many wildlife population dynamic and management 
models include the effects of climatic variables (e.g., Dennis & Otten, 2000; Colchero etal., 2009) 

3,4" 



• "The use of these climatic variables: the average annual precipitation and winter 
temperature. Without any justification, authors could have also included others variables in 
the models, as climatic variables do not account for the pattern described, but this doesn't 
necessarily mean that climate hasn't got any effect on the guanaco population growth rate. 
Moreover authors should explain if these two climatic variables are selected as indicators of 
primary production and harsh weather conditions. " 

To explain and justify the use of these clinnatic variables the paragraph below was 

incorporated in the section "Clinnatic variables": 
"These climatic variables were selected because some studies suggested that they have influence 
on guanaco demographic parameters: Rey et al. 2012^^ observed that after a severe drought the 
proportion guanaco yearlings/female decreased significantly, and Sarno et al. 1999'^^ detected a 
negative effect of winter snowfall on guanaco juvenile survival; because we had insufficient 
snowfall data we used the winter temperature as a proxy for winter snowfall. " 

• "As authors stated in introduction, guanaco populations have undergone overhunting 
(poaching) and competition with sheep, these being two of the major forces of their dramatic 
decline in the last century. For that reason, it is quite surprising that no variables related to 
livestock were included in the modelling procedure. The population has increased a lot in 
the studied period (x8) and just climate and/or density-dependent process do not explain 
alone such increase in my view. Therefore authors should consider including livestock 
effects, its influence on resource availability along the studied period and the mentioned 
overhunting, as this could really explain the detected trend. " 

In relation with the "poaching" the following paragraph was incorporated in section "Study 
area" to justify why this variable was not considered in our study: 

"In contrast to the continent, the puma (Puma concolor), the main predator of guanacos, is 
absent in the island; since 1977 the guanaco hunting has been controlled by Chilean 
National Forest Corporation (CONAF, according to its Spanish acronym) and the Chilean 
Agricultural and Livestock Service (SAG, according to its Spanish acronym)." 

In relation to the suggestion to incorporate livestock in the analysis, we re-analyzed our 
multiple regression incorporating the variable livestock as sheep population. We 
re-processed the regressions incorporating the sheep population as a new independent 
variable as a possible predictor of the guanaco population growth rate (I); so our new 
analysis increased from 5 to 6 independent variables: sheep population, guanaco 
population size, winter temperature, winter temperature with a lag of 1 year, annual 
precipitation, annual precipitation with a lag of 3 years. Our main results didn't change and it 
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was confirmed that only the guanaco population size seems to have an effect on the 
population growth rate (I). These changes can be found in the "Methods" and in the 
"Results" sections. 

• "The guanaco sampling method has changed during the studied period (different band 
widths). Besides, were the transects chosen based on the topography and vegetation cover 
of the ranch to ensure that aii areas were visibie to census-talkers aiong their transects? I am 
not sure if that couid be a source of error, but authors shouid discuss this point. " 

In the Section "Guanaco population sampling" we added the following paragraphs where we 
explain the motive of the change and justify the transects chosen: 

"The sampling methodology was changed by an expert recommendation made in 1995 as in 
year 2000 the fixed and the undefined band width methods were applied simultaneously, we tested 
both results and found no significant statistical difference between them. " 
"Despite randomly selected transects are recommended^"^, preexisting roads were used because 
according to Soto (2010)^^ the existing system of roads is an adequate sample of the area." 



''Moreover, I recommend authors to re-analyze their data using the DISTANCE program 
(Buckland et al., 2004) to estimate population size over the period, as the analyses 
performed by this package seems much more robust that the one employed by authors. If 
they do so, authors could also include topographic variables in the models to see whether 
the resulting models change. The reason is that the authors do not explain ifguanacos 
occur in both biomes (steppe and forest), because if that was true, detection ofguanacos 
was completely different in both ecosystems. Therefore the band width of 500 or 1000m is 
not valid for forests. " 



In the Section "Guanaco population sampling" the following paragraphs was incorporated to 
justify the methodology used: 

"Raedeke ( 1978) claims that the Leopold method is valid when the following conditions 
are satisfied: (1)the road systems must be an adequate sampling of the study area, (2) the 
network road must be randomly distributed, and (3) the animals included in the sampling 
must be randomly distributed in relation to the observer's route; and he considers that in the 
south of the Tierra del Fuego Island these assumptions are fulfilled. Soto 2010^^ compared 
population estimations by the Leopold and Distance methods and found that the value of 
the means estimate by Leopold methods fall within the confidence intervals estimated by 
Distance; additionally, as all sampling periods used the same methodology whatever bias 
might exist in the estimation of the mean abundance using the Leopold method, the relative 
changes among years (and thus the temporal trend) will be adequately represented; thus, 
we conclude that the Leopold's method seems adequate for our purposes " 



In relation to the forest and the sampling methodology, as the area occupied by forests is 
low and it has some degree of intervention, the visibility of the guanaco is adequate, hence 
we considered that the forest is a minor component and that does not affect the efficiency of 
the sampling procedure. In the section "Study area" the next paragraph was incorporated: 
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"The forested area is about 8.8% of the total study area, with different degrees of forest 
clearance that offers adequate visibility forguanaco sampling." 

"Regarding the population size estimated along the 36 years (Supplementary table), 
important (huge) differences are found, for example between 2004 and 2005 or between 
2006 and 2007, or between 2007 and 2008. My concern is if such "jumps" could be related 
to the census surveyed. Please discuss this point." 

In the "Discussion" section we added tine following discussion as required by the reviewer: 
"The results of the guanaco population sampling suggest a certain trend in the population 
size, with a sort of exponential growth in the first few years, becoming more variable as the 
population grows, with more marked fluctuations during the last 10-12 years; this may imply 
that the population is becoming stabilized, possibly approaching its carrying capacity. On 
the other hand, the sampling results of those last 10- 12 years show abrupt "jumps" in some 
years (for example 2004-2005) that may be the consequence of the mobility ofguanacos 
from the Cameron ranch to neighboring sites and vice versa; since the transect sampling 
does not identify individuals the effects of possible local displacements could not be 
considered." 



• "The authors used a stepwise regression when other approaches such as logistic 
regressions, Ricker or Beverton-Holt models could give more suitable results. Please 
explain or justify your decision. " 

Our goal was to carry out a preliminary analysis that would guide us in the study of the 
population regulation process in the guanaco. As there are studies in ungulate that have 
used simple multiple regression analysis and obtained good/satisfactory results, we 
considered that it was of interest to test this technique as a first step in such an analysis. 
In the manuscript the following paragraph was incorporated in the "Discussion" section to 
justify our decision: 

"Contrary to what was expected based on the bibliography of ungulate's population 
dynamics, weather variables do not seem to influence the density-dependent population 
growth rate process of the guanaco population of the Cameron ranch. Our aim was to carry 
out a preliminary analysis that would help identify some of the regulation processes in 
guanaco population; so we used a simple multiple regression analysis as a first step in that 
direction; nevertheless, our conclusion is that in order to test the effects of climatic variables 
on population regulation of large ungulates as the guanaco, the use of a population dynamic 
model would be recommended, for they seem to account better for the interaction between 
density-dependent processes and weather variables than a simple regression between the 
latter and guanaco population growth rate. " 

• "The selection of the best model among the candidates should be done byAkaike's 
information" 

This suggestion was considered and the respective Akaike's values were calculated. 
Minor comments: 
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1 . "Refer to climatic variabies or ciimatic effects, instead of ciimatic covariabies. " 



The recommendation was applied. 

2. "Expiain more in detaii Corani & Gatto's studies. " 

With the previous changes we realize that the citation of Corani and Gatto (2007) is not 
necessary any more. Mainly because the essence of that paper is the proposal of a 
methodology that is not strictly relevant to our work. The phrase with this citation was 
deleted. 

3. "Piease remove information reiated to metaboiic studies between guanacos and stieep ( 1st 
paragrapti of study area section)." 

The suggestion was carried out. 

4. "Piease add Nottiofagus spp among ttie bracl<ets after deciduous forest (study area section) 
if it is tiie case." 

This information is already incorporated. 

5. "Give more information about iivestoci< ranctiing, density, and size and tiie relation of 
iivestoci< rangers witti guanacos along the studied period (study area section). " 

This information was already incorporate in the section "Sheep population" 

6. "Are the guanacos found both in the steppe and the forest biomes?" 

Yes, the guanacos use both biomes, but they mostly use the open habitat (i.e. steppes and 
meadows). 

7. "Please re-write the paragraph related to guanaco sampling through transect method in the 
consecutive periods (it is not clear enough). " 

The paragraph was improved. 

8. "Please remove information related to binoculars and GPS brands (did the authors use the 
same gadgets across 36 years?). " 

The recommendation was carried out. 

9. "Explain the interest of using time lags (T) as covariabies in the model." 

The next paragraph was incorporated in section "Population analysis": 

"The inclusion of time lags is convenient because it may detect possible cumulative climatic 
effects on ungulate population growth rates as it was observed in, e.g., deer and moose 
and in vicunas^ '^". 
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1 0. "Please add the Statistica version employed. " 

The Statistica tools was replaced by language R (Version 0.97.449 - ©2009-2012 RStudio, 
Inc.) to obtain the value of Akaike and satisfy the item 3.2 (see above). 

1 1 . "Please move explanations of the analyses perform from results to material & methods 
section." 

The recommendation was carried out. 
Competing Interests: No competing interests 




Ricardo Baldi 

Terrestrial Ecology Research Unit, Centre Nacional Patagonico, Puerto Madyrn, Argentina 



Approved with reservations: 20 November 2013 

Referee Report: 20 November 2013 
doi:1 0.5256/f 1 000research.2588.r2048 

The authors hypothesized that population density and climatic variables affect the growth rate of the 
population of guanacos in the study area. However, as it is detailed under the "Study area" subsection, 
the domestic sheep share the area with the guanaco and substantial changes in stocking rate have been 
reported for the ranch during the study period. It is known that sheep and guanacos select similar diets 
across different parts of the range they share, and usually sheep presence and abundance explain a 
substantial amount of variation in the abundance of guanacos, included the study area comprising the 
Estancia Cameron as described by Raedeke (1979). In his work, Raedeke described that sheep are 
managed "rotationally" around the property, taken seasonally to different areas in winter and summer, and 
this makes guanacos move seasonally to the forest habitat. Therefore, it would be expected that 
variations in sheep abundance, not only seasonally but in total numbers through time, could affect the 
estimates of guanaco abundance. In my opinion, this is a major issue in the design of this work. Also, as 
explained below, the survey design of this work could play a key role in the quality of the estimates 
obtained, which is reflected also in the interpretation of results. All of these do not mean that density 
dependent processes are not likely to occur, but to my understanding it is not possible to rule out other 
processes to be conclusive. 

Guanaco population sampling 

The surveys were conducted at a 2,000 km^ ranch where the habitat is a mosaic of steppe and forest 
biomes. The authors indicate that they surveyed guanacos using the "transect method with a variable 
width" from 1 977 to 2000. Is this the distance sampling? If so, it is based on measuring the perpendicular 
distance from the transect to the object. Also, it is says that the surveying methodology was changed to 
"fixed width band" in 2001 and subsequent years. Some comments and questions on these issues 
regarding survey design are: 

1 . Were both the steppe and forest habitats surveyed using the same methodology? Transect lines to 
observe individuals and make direct counts are commonly used in open habitats, but not in the 
forest. However if this was not the case, the authors need to account for the variation in the 
probability of detection as it is expected to vary significantly between such contrasting habitat 



Page 16 of 21 



FlOOOResearch 



FIOOOResearch 2014, 2:210 Last updated: 21 AUG 2014 



types. 

2. There was a change in the surveying method, which also brings the point of detection probability. 
While the line transect method assumes that the probability of detecting objects on the line 
(distance=0) is p=1 , and as objects are far from the line the probability of detection follows a 
certain declining function, the main assumption of the strip transect method is that the probability of 
detecting objects within the strip is 1 (i.e. no objects are left undetected). This variation in the 
probability of detection according to the methodology applied during each period must be 
discussed. Also, the ways that the density estimates were obtained are unclear and must be 
specified (i.e. distance sampling, other calculation methods). 

3. Regarding the estimate of population size, several considerations can be made about the 
extrapolation of absolute densities to population numbers. First, the spatial pattern of the transects 
surveyed could be determinant of biased estimates, if the objects are not randomly distributed in 
relation to the survey line (i.e. if the transect follows a geographic feature in such a way the animals 
are either attracted or deterred). Perhaps the transects are located across certain areas - for 
example the steppe - but not the forest and thus the estimate could be biased towards the density 
calculated for the surveyed portion of the total area. This is especially important at the time of 
extrapolate densities to obtain population size numbers, as densities are strictly valid for the areas 
effectively surveyed while traveling along the transect lines. In any case, it would be useful to 
provide the distribution of surveyed lines across the area. The second point I find relevant to 
discuss regarding survey design is the proportion of the area effectively surveyed. How many km^ 
were comprised by the lines and strips surveyed? Again, the extrapolation of densities to 
population size at a larger scale must account for possible limitations related to the survey design 
and effort invested. 

Results and discussion 

Major fluctuations of the guanaco population occurred during the last 10-12 years analyzed, as can be 
found in the additional data file provided, and also pointed out by the authors. But this variation could well 
be due to the changes in methodology at the time of the population surveys. As said above, it is 
necessary to provide detailed information on the type of survey conducted, the differences with the 
previous method and the way the density estimates were obtained. It would be informative if the authors 
provide the errors associated with the population estimates. 

As proposed above, sheep densities could be affecting guanaco population trends, but the design of the 
work does not allow accounting for this factor. In my opinion, sheep densities must be included as an 
explanatory variable as there is evidence they affect guanaco numbers. Also, as the authors concluded 
that in this study the guanaco populations was near carrying capacity, and this must be common to both 
guanacos and sheep, thus sheep must be considered in the analysis as its density was markedly variable 
during the study period. Another factor that should be considered at the time of the discussion is the 
possibility of guanacos moving from the surveyed surroundings or even from Cameron to neighboring 
sites, as they may respond to changes in the spatial availability of resources by occupying different sites. 
In summary, the marked fluctuations in population size require the examination and discussion of other 
possible factors before concluding that density dependent processes would drive population trends. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard, however I have significant reservations, as outlined 
above. 
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1 Comment 

Author Response 

Maria Zubillaga, CONICET-UNLP, Argentina 
Posted: 14 IVlay 2014 

• "The authors hypothesized that population density and climatic variables affect the growth 
rate of the population ofguanacos in the study area. However, as it is detailed under the 
"Study area" subsection, the domestic sheep share the area with the guanaco and 
substantial changes in stocking rate have been reported for the ranch during the study 
period. It is known that sheep and guanacos select similar diets across different parts of the 
range they share, and usually sheep presence and abundance explain a substantial amount 
of variation in the abundance of guanacos, included the study area comprising the Estancia 
Cameron as described by Raedeke (1979). In his work, Raedeke described that sheep are 
managed "rotationally" around the property, taken seasonally to different areas in winter and 
summer, and this makes guanacos move seasonally to the forest habitat. Therefore, it 
would be expected that variations in sheep abundance, not only seasonally but in total 
numbers through time, could affect the estimates of guanaco abundance. In my opinion, this 
is a major issue in the design of this work. Also, as explained below, the survey design of 
this work could play a key role in the quality of the estimates obtained, which is reflected 
also in the interpretation of results. All of these do not mean that density dependent 
processes are not likely to occur, but to my understanding it is not possible to rule out other 
processes to be conclusive. " 

In relation to the suggestion to incorporate livestock in the analysis, we re-analyzed our 
regression incorporating the variable livestock as sheep population; we re-processed the 
regressions incorporating the sheep population as a new independent variable of the 
guanaco population growth rate (I); so our new analysis increased from 5 to 6 independent 
variables: sheep population, guanaco population size, winter temperature, winter 
temperature with a lag of 1 year, annual precipitation, annual precipitation with a lag of 3 
years. Our main results didn't change and it was confirmed that only the guanaco population 
size seems to have an effect on the population growth rate (I). These changes can be found 
in the "Methods" and in the "Results" paragraphs. 

• "The surveys were conducted at a 2, 000 km^ ranch where the habitat is a mosaic of steppe 
and forest biomes. The authors indicate that they surveyed guanacos using the "transect 
method with a variable width" from 1977 to 2000. Is this the distance sampling? If so, it is 
based on measuring the perpendicular distance from the transect to the object. " 

No, the sampling methodology did not use the program Distance. The citations of Buckland 
et al. (1993 and 2001) were removed to avoid possible misinterpretations with the sampling 
methodology used.The method used to estimate the population size was the King method 
modified by Leopold (1933) as described in Raedeke (1978). In the manuscript it was 
explained and justified this method with the following paragraph: 
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"Raedeke ( 1978) claims that the Leopold method Is valid when the following conditions are 
satisfied: (1)the road systems must be an adequate sampling of the study area, (2) the network 
road must be randomly distributed, and (3) the animals Included In the sampling must be randomly 
distributed In relation to the observer's route; and he considers that In the south of the TIerra del 
Fuego Island these assumptions are fulfilled. Soto 2010'^^ compared population estimations by the 
Leopold and Distance methods and found that the value of the means estimate by Leopold 
methods fall within the confidence Intervals estimated by Distance; additionally, as all sampling 
periods used the same methodology whatever bias might exist In the estimation of the mean 
abundance using the Leopold method, the relative changes among years (and thus the temporal 
trend) will be adequately represented; thus, we conclude that the Leopold's method seems 
adequate for our purposes The area effectively surveyed In each sampling period was around 
420 km^ (about 20%, of the area under study). " 

• "Also, It Is says that the surveying methodology was changed to "fixed width band" In 200 1 
and subsequent years. Some comments and questions on these Issues regarding survey 
design are: 

1 . Were both the steppe and forest habitats surveyed using the same methodology? Transect 
lines to observe Individuals and make direct counts are commonly used In open habitats, 
but not In the forest. However If this was not the case, the authors need to account for the 
variation In the probability of detection as It Is expected to vary significantly between such 
contrasting habitat types. " 

We corrected in the manuscript the years where the change in band width occurred: 
variable band width from 1 977 to 1 995, and with a fixed width band (with a maximum of 500 
m to each side of the transect) from 1 996 to 201 2. 

The area occupied by forest is low and it has some degree of intervention, the visibility of 
the guanaco is adequate, hence we considered that the forest is a minor component and it 
does not affect the efficiency of sampling. In the section "Study area" the following 
paragraph was incorporated: 

"The forested area Is about 8.8% of the total study area, with different degrees of forest clearance 

that offers adequate visibility for guanaco sampling." 

• "There was a change In the surveying method, which also brings the point of detection 
probability. While the line transect method assumes that the probability of detecting objects 
on the line (dlstance=0) Is p=1, and as objects are far from the line the probability of 
detection follows a certain declining function, the main assumption of the strip transect 
method Is that the probability of detecting objects within the strip Is 1 (I.e. no objects are left 
undetected). This variation In the probability of detection according to the methodology 
applied during each period must be discussed. Also, the ways that the density estimates 
were obtained are unclear and must be specified (I.e. distance sampling, other calculation 
methods). " 

We improved this section and the next paragraph was incorporate to explain this 

change. The population size was estimated by King method modified by Leopold (1933) as 

described by Raedeke (1978), which was also clarified in the manuscript. 
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"The sampling methodology was changed by an expert recommendation made In 1995 as In 
year 2000 the fixed and the undefined band width methods were applied simultaneously, we tested 
both results and found no significant statistical difference between them. " 

• "Regarding the estimate of population size, several considerations can be made about the 
extrapolation of absolute densities to population numbers. First, the spatial pattern of the 
transects surveyed could be determinant of biased estimates. If the objects are not 
randomly distributed In relation to the survey line (I.e. If the transect follows a geographic 
feature In such a way the animals are either attracted or deterred). Perhaps the transects 
are located across certain areas - for example the steppe - but not the forest and thus the 
estimate could be biased towards the density calculated for the surveyed portion of the total 
area. This Is especially Important at the time of extrapolate densities to obtain population 
size numbers, as densities are strictly valid for the areas effectively surveyed while traveling 
along the transect lines. In any case. It would be useful to provide the distribution of 
surveyed lines across the area. The second point I find relevant to discuss regarding survey 
design Is the proportion of the area effectively surveyed. How many km^ were comprised by 
the lines and strips surveyed? Again, the extrapolation of densities to population size at a 
larger scale must account for possible limitations related to the survey design and effort 
Invested. " 

Although the transects are not randomly distributed, the preexisting roads can be 
considered as an adequate sample of the area. This comment was incorporate in the 
manuscript. 

"Despite randomly selected transects are recommended^^, preexisting roads were used because 
according to Soto (2010)'^^ the existing system of roads Is an adequate sample of the area." 

• The second point was considered and the area effectively surveyed was estimated. In the 
"Guanaco population sampling" section the following paragraph was incorporated: 

"The area effectively surveyed In each sampling period was around 420 km^ (about 20%, of the 
area understudy)." 

• "Mayor fluctuations of the guanaco population occurred during the last 10- 12 years 
analyzed, as can be found In the additional data file provided, and also pointed out by the 
authors. But this variation could well be due to the changes In methodology at the time of 
the population surveys. As said above. It Is necessary to provide detailed Information on the 
type of survey conducted, the differences with the previous method and the way the density 
estimates were obtained. It would be Informative If the authors provide the errors associated 
with the population estimates. " 

The errors associated with the population estimates are given in Table S1 where the lower 
and upper limits of the 95% confidence intervals (CI) are shown. 

• "As proposed above, sheep densities could be affecting guanaco population trends, but the 
design of the work does not allow accounting for this factor. In my opinion, sheep densities 
must be Included as an explanatory variable as there Is evidence they affect guanaco 
numbers. Also, as the authors concluded that In this study the guanaco populations was 
near carrying capacity, and this must be common to both guanacos and sheep, thus sheep 
must be considered In the analysis as Its density was markedly variable during the study 
period. Another factor that should be considered at the time of the discussion Is the 
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possibility ofguanacos moving from tine surveyed surroundings or even from Cameron to 
neigiiboring sites, as ttiey may respond to ctianges in tine spatiai avaiiabiiity of resources by 
occupying different sites. In summary, the marked fluctuations in population size require the 
examination and discussion of other possible factors before concluding that density 
dependent processes would drive population trends. " 

This factor (sheep population) has been incorporated, and the corresponding analysis was 
carried out; see in the "Materials and methods" section "Sheep Population" and the section 
"Results". In the "Discussion" section the following paragraph was included where we 
recognized the possibility of guanaco moving as possible source of error in the sampling: 
"The results of the guanaco population sampling suggest a certain trend in the population size, with 
a sort of exponential growth in the first few years, becoming more variable as the population grows, 
with more marked fluctuations during the last 10- 12 years; this may imply that the population is 
becoming stabilized, possibly approaching its carrying capacity. On the other hand, the sampling 
results of those last 10- 12 years show abrupt "jumps" in some years (for example 2004-2005) that 
maybe the consequence of the mobility ofguanacos from the Cameron ranch to neighboring sites 
and vice versa; since the transect sampling does not identify individuals the effects of possible 
local displacements could not be considered. " 
Competing Interests: None 
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